Interaction corrections: temperature and parallel field dependencies of the Lorentz 

number in two-dimensional disordered metals. 
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The electron-electron interaction corrections to the transport coefficients are calculated for a 
two-dimensional disordered metal in a parallel magnetic field via the quantum kinetic equation 
approach. For the thermal transport, three regimes (diffusive, quasiballistic and truly ballistic) 
can be identified as the temperature increases. For the diffusive and quasiballistic regimes, the 
Lorentz number dependence on the temperature and on the magnetic field is studied. The electron- 
electron interactions induce deviations from the Wiedemann-Franz law, whose sign depend on the 
temperature: at low temperatures the long-range part of the Coulomb interaction gives a positive 
correction, while at higher temperature the inelastic collisions dominate the negative correction. 
By applying a parallel field, the Lorentz number becomes a non-monotonic function of field and 
temperature for all values of the Fermi-liquid interaction parameter in the diffusive regime, while in 
the quasiballistic case this is true only sufficiently far from the Stoner instability. 

PACS numbers: 71.10.Ay, 72.15.Gd, 72.15.Eb 



I. INTRODUCTION 

A standard result of the Drude-like theory of transport 
in disordered metals is the Wiedemann-Franz law 1 relat- 
ing the (Drude) thermal (k d ) and electrical (o~ D ) conduc- 
tivities via the Lorentz number Lq: 
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where T is the temperature in energy units (kg = 1) and 
e is the electronic charge. "Drude-like theory" means 
that two assumptions are made in order to calculate 
the transport coefficients: 1. the electrons do not in- 
teract with each other; 2. the scattering of the elec- 
trons onto the impurities is elastic. 2,3 While it was shown 
long ago 4 that the interplay of electron-electron inter- 
actions and disorder leads to logarithmically divergent, 
temperature-dependent corrections to the electrical con- 
ductivity at low temperatures T -C h/r (r is the mean 
free time for the impurity scattering), it is only recently 
that such effects have been correctly evaluated at higher 
temperatures 5 and for the thermal transport. 6-8 In par- 
ticular early calculations 9,10 of the interaction corrections 
to the thermal conductivity arrived at contradictory re- 
sults, due to technical difficulties in the proper construc- 
tion of the energy current density operator (both in the 
diagrammatic technique and in the kinetic equation ap- 
proach). This issue has been resolved in Ref. 6, where the 
local form of the collision integral for the kinetic equation 
is also presented. 

In deriving the quantum kinetic equation, a proper de- 
scription of the disordered Fermi-liquid is obtained by in- 
troducing bosonic soft modes (interacting electron-hole 
pairs) which contribute to the energy transport but, be- 
ing neutral, not to the charge transport. For interac- 
tion in the triplet channel these bosons have a total spin 
L = 1; this spin degree of freedom is affected by the mag- 
netic field: the description of such effects is a central part 



of the present work. By extending the results of Ref. 6, I 
analyze in detail, for two-dimensional systems, the tem- 
perature and parallel magnetic field H dependencies of 
the "generalized" Lorentz number L, defined as 



L(T, H) = k(T, H) /a(T, H)T , 



(1.2) 



where, due to the electron-electron interaction correc- 
tions, both conductivities are temperature- and field- 
dependent (a similar analysis for zero-dimensional sys- 
tems - open quantum dots - is presented in Ref. 11; the 
parallel field dependence of a is considered in Ref. 12). 
Because of difficulties in accurately measuring the elec- 
tronic thermal conductivity, very few experiments have 
been performed in two-dimensional systems with regards 
to the thermal transport - for example the Wiedemann- 
Franz law was found to hold 13 in a Si MOSFET sample 
within the experimental accuracy, and the validity of this 
law was checked for the weak localization correction. 1 
One of the difficulties in determining k is in separat- 
ing the electronic contribution to the total thermal con- 
ductivity from the phonons' contribution; however, by 
measuring the thermal conductivity in the presence of 
a magnetic field it may be possible to extract the elec- 
tronic field-dependent part, as done e.g. for cuprate 
superconductors. 15 Since new methods for measuring the 
thermal conductivity in thin films are being explored, 16 
the study of the field dependence of L could be experi- 
mentally relevant. 

The paper is organized as follows: in the next section 
I examine the temperature dependence of the Lorentz 
number L to identify different regimes as a function of 
the dimensionless parameter Tt/Ti and to discuss the 
various approximations involved. In Sec. Ill I present 
the results for the dependence of L on the parallel mag- 
netic field. The derivation of these results is given in 
Sec. IV. After the conclusions, I briefly consider in Ap- 
pendix A the field-dependent correction to the specific 
heat. Appendix B contains some mathematical details, 
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and Appendix C a discussion of the electrical magneto- 
conductivity in parallel field. 



II. TEMPERATURE DEPENDENCE OF THE 
LORENTZ NUMBER. 

This section contains the result for the temperature 
dependence of the Lorentz number L. It is convenient 
to separate L into a "Wiedemann-Franz law" part Lo 
[Eq. (1.1)] and a "violation" part SL as follows: 



L(T) = L + Sh(T) , 



(2.1) 



and for the correction SL to distinguish the contributions 
due to interactions in the singlet and in the triplet chan- 
nels: 



SL = SL S + SL f 



(2.2) 



For clarity, the two terms are considered separately. The 
results given below are derived in Sec. IV A within the 
quantum kinetic equation approach - this is a perturba- 
tive approach with 1/g as the small parameter, where 
g = a/(2e 2 /h) >• 1 is the dimensionless conductance; it 
assumes the validity of the Fermi-liquid picture, which in 
turns requires T <§C Ep, with Ep the Fermi energy. 



A. Singlet channel. 

The singlet correction SL S is, with logarithmic accu- 
racy: 



5L S _ 1 
Lo TTg 



gi (2irTT/h) In (r s E F /T^ 
- i 92 (ttTt/Ti) In (l + (ft/2^Tr) 2 ) (2.3) 

where the functions g\ and gi, given in Eqs. (4.4), de- 
scribe the crossover form the low temperature diffusive 
regime to the higher temperature quasiballistic one, and 
r s is the "gas parameter" characterizing the interaction 
strength: 
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(2.4) 



with vp the Fermi velocity and e the dielectric constant. 

In the diffusive regime T <g; H/2ttt 7 both gi and gi 
tend to 1 and Eq. (2.3) reduces to: 
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Lo 27rg V 2T 



(2.5) 



where D = tv f /2 is the diffusion constant and k — 
2-Kve /e is the 2D inverse screening radius (where v = 
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FIG. 1: Relative correction to the Lorentz number as a func- 
tion of temperature in the diffusive regime. Solid lines are 
given by Eq. (2.3), while dashed lines by the approximate 
expression (2.5). For thick lines r 3 — 0.1 and for thin ones 
r s = 1; the conductances are for each pair of curves (top to 
bottom): g = 100, g = 400 and g = 1000. 



m/ir is the 2D density of states). Compared to the 
Altshulcr-Aronov interaction correction to the electrical 
conductivity, 4 this logarithmic correction to the Lorentz 
number has a completely different physical origin: it 
arises from the energy transported over long distances 
by the neutral bosonic soft modes of the interacting elec- 
tron system. 6 At low temperatures, this additional chan- 
nel for the energy transport (as compared to the charge 
transport) leads to an increase in the thermal conductiv- 
ity over the electrical conductivity and therefore to an 
enhancement of the "generalized" Lorentz number (1.2). 

In Fig. 1 1 plot the relative change of the Lorentz num- 
ber 5L S /Lo, Eq. (2.3), as a function of 2hTt/Ti for three 
conductances (g = 100, 400 and 1000) and for two val- 
ues of the interaction strength (r s = 0.1 and r s = 1); 
for comparison the curves obtained using the approxi- 
mate expression (2.5) are also shown. From the figure 
and the dependence on r s in Eq. (2.3) it follows that the 
deviation from the Wiedemann-Franz law grows with the 
interaction strength. 20 For low conductances and temper- 
atures the (positive) change in the Lorentz number is of 
the order of a few percent; unfortunately the uncertainty 
in measurements of the thermal conductivity in metallic 
films 16 is also of this magnitude, making a comparison 
with the present theory pointless. 

As the temperature increases, the inelastic collisions 
between the electrons and the bosons tend to inhibit 
the energy transport more efficiently. The quasiballistic 
regime is reached in the temperature range 



h/2irT «T« T° b , 



(2.6) 
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plotted (dashed lines). In agreement with the above dis- 
cussion, comparison of the dashed and solid lines shows 
that Eq. (2.9) deviates significantly from the full expres- 
sion (2.3) at low conductance, while there is good agree- 
ment at high conductance. The intersections between 
the solid lines and the upper thin dotted curve are at 
T = 0.1T* b ; it is evident that the region of validity of 
the quasiballistic approximation grows with the conduc- 
tance. For all conductances the (negative) correction can 
be a few percent; in Ref. 13 the Lorentz number was 
measured 21 in a 2DEG and found to be slightly smaller 
than Lo, in qualitative agreement with the predictions in 
the present work. However the uncertainties are of the 
same order of the calculated effect and hence no quanti- 
tative comparison is possible. 



FIG. 2: Relative correction to the Lorentz number as a func- 
tion of temperature in the quasiballistic regime. Solid lines 
are given by Eq. (2.3) with r 3 = 0.1, while dashed lines by 
the approximate expression (2.9). The conductances are (left 
to right): g = 100, 400, 1000, 1600 and 2500. Thin dotted 
lines are calculated so that they intersect the solid lines at 
T = 0.1T q s ft (top curve) and T = 0.2T q s fc (bottom curve). 



B. Triplet channel. 

For the triplet channel interaction correction, I con- 
sider separately, for simplicity, the diffusive and quasi- 
ballistic regimes; in the former case the correction is: 



where T^ b is the solution of: 
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for large conductances this gives: 
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(2.7) 



(2.8) 



In this regime, the dominant contribution to the singlet 
correction (2.3) is: 



SL. 
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According to condition (2.6), this expression is applica- 
ble if Ty b ^> Ti/2-kt; this can be satisfied only for large 
enough conductances. For example at g — 14 I find (nu- 
merically) T° b w 0.UE F w 10h/2TTT, and at g = 720, 
Ty b w O.OlEp w 50h/2TTT; in the latter case, and for 
larger conductances, Eq. (2.9) can be expected to have 
a sufficiently large range of validity, while in the former 
(and generally for small conductances) there is no quasi- 
ballistic regime. At temperatures of the order of T^ b the 
energy transport becomes truly ballistic in nature, as the 
dominant processes responsible for the relaxation of the 
energy current are the inelastic electron-boson collisions 
and not the electron-impurity collisions. Although the 
high temperature regime T > T^ b is not considered here, 
it can be treated within the kinetic equation approach. 18 
Fig. 2 shows the relative change of the Lorentz number 
<5L S /L as a function of 2ttTt/H for five different conduc- 
tances (solid lines); for the lowest and highest conduc- 
tances considered, the approximate result (2.9) is also 



■ln(l + K) 
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(2-10) 

where the first term on the right hand side is again due 
to the bosonic energy transport, and the second one 
originates from the interaction-induced energy depen- 
dence of the elastic cross section. While the sign of this 
temperature-independent correction is determined by the 
sign of the Landau Fermi- liquid constant Ffi , its contri- 
bution to the total correction SL [Eq. (2.2)] is generally 
small 22 and the overall positive sign of SL at low enough 
temperatures is determined by 5L s d . I do not plot sep- 
arately the contribution (2.10) to SL, since its effect is 
simply to shift upwards (downwards) the curves in Fig. 1 
for Fq > {F$ < 0). 

In the quasiballistic regime the correction reads: 



3SL qb ( 
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(2.11) 



with 5Lq b given in Eq. (2.9). As for the singlet chan- 
nel, this negative correction originates from the inelastic 
electron-boson collision, and similarly to the singlet chan- 
nel correction, the validity of this expression is limited at 
high temperatures by T l qb - the equation defining this 
quantity is obtained by multiplying the left hand side of 
Eq. (2.7) by 3(F CT /1 + F CT ) 2 . It is evident that, when 
the quasiballistic approximation is applicable, plotting 
the sum of Eq. (2.11) and (2.9) would give Fig. 2 with a 
rescaled vertical axis, as for any value of F£ the triplet 
channel contribution enhances the singlet channel correc- 
tion. 
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III. PARALLEL FIELD DEPENDENCE OF THE 
LORENTZ NUMBER. 

In this section I present the results for the parallel mag- 
netic field dependence of the Lorentz number. The par- 
allel field H affects the electrons by shifting the energy 
levels by the Zeeman energy 



Ez = 9l^bH . 



(3.1) 



where is the Lande g-factor and fis the Bohr mag- 
neton. The Lorentz number depends on H only through 
this energy and the renormalized Zeeman energy E* z : 



E*= Ez 
z 1 + F£ 



(3.2) 



As it is the case for other transport properties (e.g. 
the magneto-conductivity), it is convenient to consider 
the deviation AL of the Lorentz number from its zero- 
field value: 



AL(T, H) = L(T, H) - L(T, 0) . 



(3.3) 



Once again I address separately the diffusive and quasi- 
ballistic regimes; in both regimes the system is assumed 
to be far from the full polarization, i.e. Ez <C Ep. The 
derivation of the results can be found in Sec. IV B. 



A. Diffusive regime. 

For T <C h/2nT, the field-induced change in the 
Lorentz number is: 



AL„ 
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(3.4) 



with the dimcnsionless functions l\ and I2 defined in 
Eqs. (4.16)-(4.17). Similarly to Eq. (2.10), the terms with 
the numerical prefactor 3/2 are due to the bosonic en- 
ergy transport, while the remaining ones, proportional to 
1/Fq, originate from the energy dependence of the elastic 
cross section. The structure of this expression as the dif- 
ference of terms which depend on different energy scales 
(i.e. Ez and E* z ) can be traced back to the structure 
of the quantum kinetic equation, in which the bosonic 
contributions to the collision integral always appear as 
differences between a soft mode part and a "ghost" part. 6 
In the weak field limit E Z ,E* < 2nT Eq. (3.4) be- 
comes: 



AL d 



7Tg 



Ez 
2^f 



with 



fd(x) = 



x(4 + 3x) 
(1 + z) 2 



(3.5) 



(3.6) 
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FIG. 3: Relative deviation wgALtj/Lo of the Lorentz number 
from its zero-field value in the diffusive regime. Solid lines 
are given by Eq. (3.4) with, from top to bottom: Fq — —0.7, 
-0.4, -0.2, 0.2 and 0.4. For Fq = -0.7 and 0.4 the approxi- 
mate expressions (3.5) (dashed lines) and (3.7) (dot-dashed) 
are also shown for comparison. 



In the opposite case Ez,E z 3> 2irT the approximate 
formula is: 



AL ri - 



-— ln(l 

TTg 



Lo 

or equivalently: 



F n 



2 

3^g 



1-— ln(l + ^) 
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(3.7) 



(3.8) 



with SL d given in Eq. (2.10). This result can be explained 
as follows: 12 in the diffusive regime the correction is dom- 
inated by processes with small energy and momentum 
exchange, and in the strong field the bosonic modes with 
non-zero spin projection become gapped with the gap 
given by the Zeeman energy. Therefore the contributions 
due to these modes must drop out from the total correc- 
tion to the Lorentz number: this is why the correction 
(3.7) partially cancels the one given in Eq. (2.10), with 
the surviving contribution originating from the modes 
with zero total spin projection. 

In Fig. 3 the relative deviation AL^/Lo, multiplied 
by irg, is plotted as a function of Ez/2irT for differ- 
ent values of the parameter Fq . At fields such that the 
Zeeman energy is larger than temperature the deviation 
becomes quickly field- independent, but near Ez ~ 2irT 
all the curves are non-monotonic; the presence of peaks 
is due to the above discussed dependence on the two 
different energies Ez and E z , and through the latter 
(and the 1/Fq prefactors) the peaks' positions depend 
on Fq . The temperature dependence of the deviation at 
fixed field can also be read from this graph by following 
the curves from the right (low temperature) to the left 
(high temperature): the deviation is temperature inde- 
pendent at low temperatures T <C Ez/2ir and displays 
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a power-law decay (~ T~ 2 ) at high temperatures; again 
the non-monotonic behavior characterizes the intermedi- 
ate regime. 



B. Quasiballistic regime. 

Here I consider the quasiballistic regime Ti/2-kt «T« 
T* 6 , with T* 6 defined after Eq. (2.11). In this case I find: 



AL 



qb 
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(3.9) 

with 7 3 given in Eq. (4.26). For E* < 2nT the result 
takes the form: 
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while with logarithmic accuracy the large field limit 
E* > 2ttT is: 



AL 



qb 
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2 /27tTt\ : 
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(3.12) 



Comparison of Eq. (3.12) with Eq. (2.11) shows that the 
partial cancelation that was found in the diffusive limit is 
also realized in the quasiballistic one, but with an impor- 
tant difference: the gapped modes still contribute to the 
total correction to the Lorentz number because the qua- 
siballistic corrections are dominated by the inelastic scat- 
tering with large momentum exchange. The gap there- 
fore excludes the low energy (E < Ez) contributions, but 
the higher energy ones {Ez < E < Ef) are still relevant. 

Fig. 4 shows the field dependence of AL q b, Eq. (3.9), 
by plotting 

AL 9 fc / 3 (toTrV 
U / tog\ h J 

as a function of Ez /2irT; for comparison the approximate 
result (3.12) is also plotted. 23 In the inset the low-field 
behavior of Eq. (3.9) is compared to Eq. (3.10); except for 



the case FX 



-0.7, all curves are non-monotonic. The 



threshold value F° h above which AL^ is a non-monotonic 
function can be found by requiring fqb(Ff h ) = 0; this 
gives F° h ~ —0.679. Although the present results are 
not valid close to the Stoner instability (see footnote 22 ), 
they suggests that as — > — 1 the relationship between 
energy and charge transport properties can be qualita- 
tively altered compared to the weakly interacting case. 
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FIG. 4: Relative deviation of the Lorentz number from its 
zero-field value as a function of the Zeeman energy in the 
quasiballistic regime (normalized as explained in the text) 
for different values of the Fermi-liquid parameter, as labeled 
(the curve partially covered by the inset corresponds to Fq = 
—0.7). Solid lines are given by Eq. (3.9), while dashed lines 
by the approximate expression (3.12). The inset shows details 
of the low-field regime, with the dashed curves corresponding 
to Eq. (3.10). 



0.03 



0.02 



AL, 



'qb 



Lo 



0.01 



0.00 



-0.01 




0.0 0.5 



1.0 1.5 2.0 

2nT/E z 



FIG. 5: Relative deviation of the Lorentz number from its 
zero-field value as a function of the temperature in the qua- 
siballistic regime (normalized as explained in the text) for 
different values of the Fermi-liquid parameter. Proceeding 
clockwise near the origin, starting with the steepest curve, 
the parameters are: = -0.7, -0.5, -0.3, -0.2, 0.4 and 
0.2 



For completeness, I consider in Fig. 5 the temperature 
dependence of AL q b by plotting 



AL 



qb 



1 3 
2^g 



E z t 



as a function of 2irT/Ez in the low to intermediate tem- 
perature regime. 
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IV. DERIVATION. 

This section is devoted to the calculation of the in- 
teraction corrections to the Lorentz number using the 
formalism of Rcf. 6. In the absence of the magnetic field, 
one can use directly the results of that reference, while 
a generalization is needed for the parallel field 
discussed in Sec. IV B. From now on, I set h = 1. 



A. Temperature dependence. 

The results presented in Sec. II are a straightforward 
consequence of the findings of Ref. 6, where it is shown 
that for two-dimensional systems the thermal conductiv- 
ity can be written as: 



+ Ak. 



Here 



= L ctT 



(4.1) 



(4.2) 



follows the Wiedemann-Franz law with Lo defined in 
Eq. (1.1) and the electrical conductivity a includes the 
interaction corrections. The additional term Ak = Ak s + 
3 Ak* is given by the sum of the singlet and triplet channel 
corrections. The former was calculated with logarithmic 
accuracy: 24 



the asymptotic forms of the functions g\ and g 2 given af- 
ter Eqs. (4.4). The condition (2.7) is obtained by equat- 
ing the (absolute value of the) correction (2.9) to the 
non-interacting Lorentz number L [Eq. (1.1)]; the results 
are not valid at temperatures higher than T^ b because in 
solving the kinetic equation it was assumed that the im- 
purity scattering is the dominant process contributing to 
the energy relaxation rate - see the discussion at the end 
of Sec. 6.2 of Rcf. 6; however, the kinetic equation itself 
is still valid. 

For the triplet channel, the correction Ak* was calcu- 
lated in the diffusive and quasiballistic regimes: 



A 4 = -777 



T 
18 



1 



^0 



^ln(l + K) 



for Tt < l/2?r, and 
2tt 2 



A< 6 



15 



T {Tt) In 



~T~ 



FR 



1 



(4.6) 



(4.7) 



for Tt > 1/2tt. Multiplying Eqs. (4.6)-(4.7) by 3 and 
using Eq. (4.5) and the definition (2.1) of SL, the results 
(2.10) and (2.11) are obtained; the factor of 3 arises from 
the summation over the three projections of the total 
spin, which in the absence of magnetic field contribute 
equally to the thermal conductivity. 



B. Parallel field dependence. 



Ak s 



24 
2tt 2 



g 2 (ttTt) In 1 + 



(27tTt) 2 



(4.3) 



_T(Tr) 2 In 
15 1 ; \ T 



with the functions g\ and g 2 given by: 
3 f 1 



9i{x) 



2tp'(l/x) -x 2 



(4.4a) 



where tp' is the derivative of the digamma function, and 

g 2 (x) = ^x 2 + ^ gi (x)-^. (4.4b) 

Note that g x (0) = g 2 (0) = 1, and g x {x) ~ 3/x for x < 1. 

Using the definitions (1.2), (4.1) and (4.2) I find at first 
order in 1/g: 



UT) 
Ln 



Ak 
a D T 



(4.5) 



Then from the definition (2.1) and Eq. (4.3) one arrives 
at Eq. (2.3). The limiting expressions given in Eqs. (2.5) 
and (2.9) are found by keeping the leading order contribu- 
tions to SL S for small and large Tt respectively by using 



To obtain the results of Sec. Ill I give here an exten- 
sion of the calculations of Ref. 6. As in Eq. (4.1), I sep- 
arate the thermal conductivity in a part which follows 
the Wiedemann-Franz law and a correction; both term 
now depend on the applied parallel magnetic field H, or 
more precisely on the Zeeman splitting, Eq. (3.1). The 
term k wf (T, H ) is straightforwardly calculated using the 
results of Rcf. 12, so one needs to consider only the cor- 
rection Ak(T,H). As discussed in e.g. Refs. 12 and 17, 
the singlet and the triplet L z = contributions to k are 
not affected by the parallel field [L z is the projection of 
the total spin along the field direction] . The effect of the 
field on the remaining L z = ±1 components of the triplet 
channel correction is to shift the frequency of the inter- 
action propagators; in the langauge of Ref. 6, the bosonic 
propagators C a (u, q: m, n 2 ; L z ) are: 

£ a (n 1 ,n 2 - 1 L z ) = fl 2 S {nTn 2 ) C (ni; L z )+ 



C (ni; L z )C (n 2 ; L z ) 



(-ic^gjy + l) C(L Z ) (4 . 8) 

C{L z )-(-i^ + ±) 1 

and the corresponding (triplet) "ghost" propagators 
C 9 (tu,q;ni 1 n 2 ; L z ) are: 25 

£ 9 (n 1 ,n 2 ; L z ) = Q 2 S (nrn 2 ) C {n 1 ;L z )+ 

\C{L Z ) (4.9) 



£ {ni;L z )£ (n 2 ;L z 



C(L Z ) 
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where 



and 



- L 2J B|) + it;, • q + l/r ^ 
C(L Z ) = y/(-i(w-L z E* z ) + l/T) 2 + (v F q) 2 . 

In the above formulas I dropped the variables co, q for 
compactness, Vi — vprii and E* z is the Zeeman energy 
renormalized by the interactions, Eq. (3.2). 

The calculation of the field-dependent thermal conduc- 
tivity proceeds now as in Ref. 6: the evaluation of the 
transport coefficients can be reduced to integrals over 
the frequency lo whose integrands consist of a distribution 
function part times a kernel part K(w); the latter is found 
after integration over the momentum q and summation 
over the total spin projections. The field-dependent ker- 
nels can be found using the expressions (4.8)- (4.9) instead 
of their zero- field counterparts; below I calculate explic- 
itly the correction Aft m to the thermal magnetoconduc- 
tivity K m {T,H) = k(T,H) - k(T,0). In other words, I 
want to write n(T,H) in the form (4.1) and define Ak„ 
as: 



An m {T, H) = Ak{T, H) - Ak(T, 0) , 



(4.11) 



AB° 



<T D 2ir 2 uj 



In 



, ,2 T?*"2 

w - E z 



(4.14b) 



As discussed above, these kernels are found by substitut- 
ing the expressions (4.8)-(4.9) for the propagators into 
the definitions of £ and B° given in Eqs. (6.9) and (6.36d) 
of Ref. 6. 26 Substituting Eqs. (4.14) into Eqs. (4.13) [and 
by a change of variable u> — ► 2itTuj] I get: 



+ 



and 



6F§ 

Ez 
2^f 



T 
T 



E% 



2irT 



E z 
2ttT 



h 



Ez_ 
2irT 



h(E) 



7T / cllO ■ 
J 



sinh ttlo 



■In 



LO 2 



(4.15a) 



where I introduced the dimcnsionlcss functions: 



(4.15b) 



(4.16) 



where Ak(T, 0) is the correction considered in the previ- 
ous section. To calculate A« m I introduce for each kernel 
K the corresponding kernel difference AK = K(H) — K(0) 
between the kernel calculated with and without the field. 
As in the preceding section, I consider separately the dif- 
fusive and quasiballistic regimes. 



and 



I 2 {E) =tt duo 



CO 



sinh 2 - 



In 



1+co/E 



l-co/E 

Eq. (4.16) can be identically rewritten as: 



(4.17) 



1. Diffusive regime. 

In the diffusive regime, the two main contributions to 
Ak come from the energy dependence of the elastic cross 
section and from the bosonic energy transport. Writing 



(4.12) 



the two terms are given by [cf. Eqs. (6.11b) and (6.36b) 
of Ref. 6]: 



Sk 



el,m 



~ 2 H Jdu)A£(uj) 



dNp 



12 duo 



and 



(4.13a) 



OO 

h(E) = 2E 2 -iJ2 n ln(l+ (S/n) 2 ) - - 



(E/nf 



+ (E/n) 2 
(4.18) 

which is useful to obtain the E <C 1 expansion, and as: 



h{E) = | In \E\ + C + -wjdu) 



sinh 2 



CO" LO 
In 1 _ — 

E 2 

(4.19) 

which gives the E » 1 asymptotic behavior; the constant 
C appearing above is: 



C = -2tt / (U- 
2 



7 + In 2ir 



sinh no) 
3 



In M 



(4.20) 



-_ C '(2)~0.99, 



r u> 3 dN P 
T dto 



with the kernels 

P 2 

A£ = - 



1 



win 



a D 2ir 2 L0 2 F °" 

+ E z \n 



to 



El 



co 2 - Ef 
to + Ez lo 



E% 



to ■ 



E, 



CO 



E* 



(4.13b) 



(4.14a) 



where 7 ~ 0.577 is the Euler constant and C'(2) ~ -0.938 
is the derivative of the zeta function evaluated at 2. As 
for Eq. (4.17), in the given form the large E limit can be 
readily obtained \h{E) = 2/3E + . . .], while to find the 
small E limit I rewrite it as: 



h{E) =7rsgn£ -IE 



+ 4^]lm 



71=1 



In 1 



E 



1 



1 - iE/n 



(4.21) 
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Note also the identity [primes indicate derivatives]: 

I[{E) +EI' 2 (E) = (4.22) 

which enables to verify that the correction 5n e i :in of 
Eq. (4.15a) vanishes in the limit Ffi — > 0, as expected. 

Using Eqs. (4.15), together with the definitions (4.1), 
(4.11) and (3.3), and dropping terms of higher order in 
1/g, one arrives at Eq. (3.4). The approximate expres- 
sions (3.5)-(3.7) follow from Eqs. (4.18)-(4.21). 

2. Quasiballistic regime. 

In the quasiballistic regime, the correction Ak is de- 
termined by the inelastic electron-boson collisions and it 
can be written as [cf. Eq. (6.39a) of Ref. 6]: 



dN P 



(4.23) 



with 
AB 1 



a D 2w 2 



T LO 



r 



1 + F% 



' p* per 



J(E;F) = -2J2 -L Z E\— { - 

L z =±l ^ { 



u 1 + F° 

(4.24a) 

{l-L z Ef 



-lnlFlF 2 



L Z E) 2 - F 2 
1 



(\-L z E) 2 -F 2 l-F 2 



(4.24b) 



Some details on the derivation of this kernel are given in 
Appendix B. Substitution of Eqs. (4.24) into Eq. (4.23) 
results in 



Ak 



m,qb 



(ilk 



-<3 



2ttT' 1 + F ct 

(4.25) 

with the dimcnsionless function I 3 defined as: 



I 3 (E;F) = n du — - 2 

sinn ituj 



H-;F 



(4.26) 



The large and small E limits of I 3 can be found by keep- 
ing the leading order terms in the expansion of the func- 
tion J(E; F). 27 In this way I obtain: 



I 3 (E;F)^--\n\E\ 



(4.27) 



for E > 1, and 



h(E;F)^ -E 2 



1 + 3F 2 , F 2 (3 + F 2 ) lnF2 



+ 



(l-F 2 ) 2 (l-F 2 ) 



(4.28) 

for E <C 1. Knowing An m ^ q b and the approximate formu- 
las (4.27)-(4.28), proceeding as in the previous subsection 
finally leads to Eq. (3.9), (3.10) and (3.12). 



V. CONCLUSIONS. 

The quantum kinetic equation approach is a powerful 
method to investigate the effects of the electron-electron 
interactions on transport in disordered metals 5 ' 6 ' 12 and 
open quantum dots. 11 Using this approach I considered 
the temperature and parallel magnetic field dependence 
of the Lorentz number in two-dimensional disordered 
metals. 

Three regimes can be distinguished as the tempera- 
ture increases: diffusive, quasiballistic and truly ballis- 
tic. In the low-temperature diffusive regime, the Lorentz 
number is enhanced above its Drude value due to the 
energy transported by neutral bosonic modes that de- 
scribe the interacting electron-hole pairs, see Eq. (2.5) 
and Fig. 1. At intermediate temperatures (the quasibal- 
listic regime) the Lorentz number is suppressed by the 
inelastic electron-boson collision, Eq. (2.9) and Fig. 2, 
with the crossover between the two regimes described by 
Eq. (2.3). The effect of the interaction in the triplet chan- 
nel is given in Eqs. (2.10) and (2.11) for the diffusive and 
quasiballistic regimes respectively. 

If a magnetic field is applied parallel to the two- 
dimensional metal, the Zeeman splitting of the electronic 
energies affects the transport properties; in particular in 
the diffusive regime the deviation of the Lorentz number 
from its zero-field value displays a non-monotonic depen- 
dence on the ratio between the Zeeman energy and the 
temperature, see Eq. (3.4) and Fig. 3. Finally, as dis- 
cussed in Sec. IIIB, in the quasiballistic regime the devi- 
ation can be either a monotonic or non-monotonic func- 
tion of both temperature and Zeeman energy depending 
on the value of the Fermi-liquid parameter. 
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APPENDIX A: PARALLEL FIELD 
DEPENDENCE OF THE SPECIFIC HEAT. 

In this appendix I calculate for completeness the cor- 
rection Sc v to the specific heat in the presence of a par- 
allel magnetic field; the general expression for 6c v is: 6 



where the energy densities u a are given by: 
u a = J dwLjb a (uj)Np(uj) . 



(Al) 



(A2) 



Here Np(uj) is the Planck distribution and b a (uj) are the 
bosonic densities of states. In the presence of the parallel 



9 



field I find: 
b a (uj;H)-b 9 (uj;H) 



Re ^ 



d 2 q 
(2tt) 2 



1 + F° \C(L 



1 



-i {u - L Z E* Z ) + 1/t 
C(L Z ) 



1 



1 



C(L z )-6 C(L z )-1/t 



with C defined in Eq. (4.10) and 

= — 



1 + itf 



(A3) 



(A4) 



Performing the integration and the summation I arrive 
at: 

(V(u,; H) - 6»(w; ff)) - (V(w; 0) - 6>; 0)) = 



8tt 2 £> 



1 + F r 



In 



o 



1 



CO 2 



In 



1 



Ef 



*2 , ,2' 



+ ttt(M - ££) ( 9(E Z - tS) - 8(E* Z 2 - to 

2 



(A5) 



o 



ttt\uj\8(E 2 z - J 2 ) 



Next, I substitute this result into Eq. (A2) and then into 
Eq. (Al); taking the temperature derivative and rescaling 
the frequency [u — > 2ttTlu] I finally obtain: 

Sc v (H) -Sc v (0) = 



1 T 

I'd 



\2irT ) 1 + F ff \2ttT J 



1 T 

7T D 



I + FZJ J3 \2TJ \ 2T 



(A6) 



^ ( 2T ) ( 2T ) (/ 2 ( 2T ) ^ 2 ( 2T 
with /i defined in Eq. (4.16) and 



fn(z) = [ 
JO 



cLu 



sinh 2 i 



(A7) 



The functions /„ can be given in terms of polylogarithms 



as: 



f 2 (z) - - + 2zlog(l-e- 22 )-Li 2 (e- 2 >) 



+z 2 (1 - cothz) 



^C(3) + 3z 2 log (1 - e- 2z ) - 3zLi 2 (e" 22 ) 



Li 3 (e- 2z ) +z 3 (l-cothz) 



(A8) 



The second line in Eq. (A6) is the correction to the 
specific heat in the diffusive limit, while the last two lines 
become dominant in the quasiballistic limit. In this limit 
and for weak interaction \Ffi\ <C 1, Eq. (A6) reproduces 
the result of Ref. 19. 



APPENDIX B: DERIVATION OF THE 
KERNEL AS 1 . 

In this appendix I briefly outline how to derive the 
kernel AB 1 given in Eq. (4.24a) starting from the results 
of Sec. 6.2 of Ref. 6. There, the exact (at linear order 
in VT) solution of the kinetic equation was given; this 
solution is unaffected by the parallel field. In the quasi- 
ballistic regime, the main corrections to the thermal con- 
ductivity were found to originate from the term in the 
bosonic distribution function defined as SN 1 and which 
can be neglected in the diffusive limit; these corrections 
are given in Eqs. (6.35c) and (6.36c) of Ref. 6. Perform- 
ing the angular integrations in those equations results in 
Eqs. (6.38a) and (6.38b) respectively; by repeating those 
calculations using the propagators in Eqs. (4.8)-(4.9) a 
similar result is obtained in the presence of the parallel 
field. The explicit expression is simply found by redefin- 
ing some of the quantities appearing on the right hand 
sides of Eqs. (6.38) as follows: for the function C, I substi- 
tute the function C(L Z ) given in Eq. (4.10); the parameter 
b' is now: 

b'{L z ) = -i(w-L z E* z ) + i; 
all the other quantities, namely b [Eq. (A4)] and 



N = VpT 



J 2 dN P Fq 
~T duj 1 + F n ff 



are unchanged. 

To arrive at the kernel AB 1 , the sum over L z need 
to be performed, along with the remaining integral over 
the magnitude of the momentum q. This integral was 
performed in Ref. 6 with logarithmic accuracy; in the 
present case however both the infrared divergence (due 
to the long-range part of the Coulomb interaction) and 
the ultraviolet divergence are absent - the latter because 
I consider here the kernel difference, the former because 
the triplet channel interaction is short-range - and I pro- 
ceed in a different manner. By rescaling all the dimen- 
sionful quantities u>, E* z , 1/t and vpq by the temperature 
T, I find that the contributions due to Eq. (6.38a) are 
smaller than those of Eq. (6.38b) by 1/TY and can there- 
fore be neglected. Similarly, to find the leading term in 
the (Laurent) expansion over 1/TY one can set 1/t — ► 
in Eq. (6.38b). Following this procedure, I get: 



TTfJ n I dp 



^Sttv 2 



xRc- 



Vp-(l-L 2 ££/a,)2 



VV - (1 - L z E* z /u;) 2 - tj^ 



1 (1 - L z E* z /wf 
P 
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where in terms of the original variables p — (vpq/oj) 2 , 
and the approximate equality indicates that I am ne- 
glecting higher order terms in 1/TV. The integral over p 
is logarithmically divergent, but the difference between 
the above expression and the similar one at zero paral- 
lel field is finite. The exact integration of this difference 
gives finally the expression for the kernel AB 1 given in 
Eqs. (4.24). 



APPENDIX C: PARALLEL FIELD 
DEPENDENCE OF THE ELECTRICAL 
CONDUCTIVITY. 

The aim of this appendix is to compare the present 
approach to the one of Ref. 12 for the calculation of the 
parallel field magneto-conductivity Act, which is given 
by: 



Ac = u D dui 



A£{u)+ AS el {u) 



d_ 

doj 



luN p {uj) . (CI) 



This formula follows from Eq. (6.8) of Ref. 6, and the 
kernels are defined in the subsequent Eq. (6.9). 

In the diffusive limit only A£ is relevant, since the ker- 
nel S el gives contributions smaller by the factor Tt< 1, 
as discussed in Ref. 6; it is straightforward to verify that 
substituting Eq. (4.14a) into Eq. (CI), the diffusive limit 
result of Ref. 12 is recovered. Viceversa, it can be shown 
that in the quasiballistic limit the kernel S can be ne- 
glected since larger corrections are due to S el - this can 
be done for example by rescaling all dimensionful quan- 
tities by the temperature, as described in Appendix B. 
Proceeding as detailed there, i.e. dropping higher order 
terms in 1/Tr, and introducing the shorthand notation: 

s(a\b) = sgn(a) - sgn(6) , 

I obtain for the kernel AS el in the quasiballistic regime: 



AS e 



AS 1 



AS 1 



(C2a) 



AS 



li 



r 7r 



<7 D 27T 2 UJ 2 



E 

L z =±\ 



s(w - L z E z \uj - L Z E* Z ) 



FZ 



x (w - L Z E* Z ) - uj s(o; - L z E z \uj) 

1 + -Tn 



(C2b) 



AS 12 ~ 



t it x - 

<j d 2tt 2 w2 ^ 

x (uj - L Z E Z ) 



=±i 



s(u> - L z E z \uj - L Z E* Z ) 



2wF% 



w(l + 2Fg) - L Z E Z 



F ^s(u-E z \uj) 



1 + F, 



(C2c) 



Eqs. (CI) and (C2) lead to the quasiballistic limit re- 
sult of Ref. 12; in particular, the sum of the first terms in 
square brackets gives the contribution denoted there by 
K2, while the remaining terms give the K\ contribution. 
Here I point out that in the low-field limit E z , E* z <C 2T, 
the quasiballistic magneto-conductance is given by: 



Act; 



where 



m = 1 



l + z 



1 1 

2(1 



1 



2z (l + 2z) 2 
-z)ln2(l + z) 



(C4) 



(l + 2z) 3 



This expression corrects the wrong definition of f(z) 
given after Eq. (14) of Ref. 12. 28 While the difference be- 
tween the two definitions is numerically small (less than 
4%) for -0.57 < Fq < 0.14, it grows rapidly outside 
this parameter range, and use of Eq. (C4) rather than its 
counterpart of Ref. 12 may be important in a comparison 
between theory and experiment. 
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The applicability of Eqs. (2.3) and (2.5) is limited to values 
of the gas parameter far from the Wigner crystal instabil- 
ity, i.e. r s < 37 - see also Sec. IIIB of Ref. 5. 
The experiment of Ref. 13 was presumably performed in 
the quasiballistic regime, since the conductivity showed a 
linear in temperature behavior (cf. Ref. 5). 
The correction (2.10) becomes large near the Stoner insta- 
bility Fq — — 1; in this regime, however, approximating 
the Landau Fermi-liquid parameter F" by its (momentum 



independent) average over the Fermi surface Ftf severely 
restricts the temperature range over which the present re- 
sults can be applied to T/E F < (1 + F ( f) 2 , see Ref. 5. 
Therefore I do not consider this regime. 
This equation gives the slope of the dashed straight lines; 
a subleading constant, calculated numerically, has been 
added to each line. 

Here an algebraic mistakes in Ref. 6 is corrected; this cor- 
rection leads to different numerical factors which, however, 
do not modify qualitatively the original results. 
The ghosts corresponding to the singlet channel are, of 
course, unchanged. 

Alternatively, one can use directly the diffusive limit ex- 
pressions for these kernels, Eqs. (6.22) and (6.24) of Ref. 6 
with the appropriately modified diffusive limit propaga- 
tors, which are obtained by shifting the frequency ui — * 
to — L Z E Z . 

The calculation of the subleading terms is more involved 
for J3 compared to the expansion of I\ or I2 and I do not 
attempt it here. 

The erroneous definition in Ref. 12 is probably the conse- 
quence of an algebraic mistake. As a further check, I have 
verified numerically that Eq. (C4) is indeed the right for- 
mula. 



